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SUMMARY 


The NTA (reference 1) Level 15.5.2/3, was used to provide 
non— linear steady— state (NLSS) and non-linear transient (NLTR) 
thermal predictions for the International Ultraviolet Explorer 
(lUE) Scientific Instrument (SI) (reference 2) . NASTRAN (references 
3 and 4) structural models (reference 5) were used as the basxs 
for the thermal models, which were produced by a straight-forward 
conversion procedure. The accuracy of this technique was sub- 
sequently demonstrated by a comparison of NTA predicts with the 
results of a thermal— vacuum test of the lUE Engineering Test 
Unit (ETU) (reference 6). Completion of these tasks was aided 
by the use of NTA pre-and-post-proc^essors and by a significant 
amount of program modification to NTA subroutines RMG and SSGHT. 


INTRODUCTION 


The NASTRAN program has quickly become a standard tool for 
the finite element structural analysis of large and/or complex 
structures. During this evolution it became apparent that the 
mathematical techniques employed in NASTRAN to solve structural 
analysis problems could be easily adapted to provide solutions 
to heat transfer problems (reference 7) , with several major 
advantages accruing, in that: 

(a) virtually the same finite element analytical model could 
be used for both cases; 

(b) with the same analytical model being used for both 
structural and thermal analyses, thermal predictions could be 
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directly applied, with no error-inducing and/or time consuming 
interpolation, to the structural model to allow the calculation 
of thermally induced deformations. 

f concept, certainly, but concepts often run 

a ou of practical details and fail to ever be used effectively 
to produce the desired result of increased productivity. The 
intent of this pager is therefore basically to increase the 
steepness of the NTA learning curve by describing in detail the 
actual techniques used in converting_an existing NASTRAN 
structural model into an executable NTA model. In support of 
tnis goal, information will be presented containing computer 
results, comparisons with test results, comments on effective 
analysis aids, and useful program coding modifications. 


THE lUE SCIENTIFIC INSTRUMENT 
Description 

The lUE SI, as shown in Figure 1, consists of a Cassegrain 
telescope section and a spectrograph section, and as a whole is 
essentially an aluminum structure approximately 3 meters long 
• diameter. The prime purpose of this instrument 

IS high resolution observation of stellar objects in the ultra- 
violet spectral region, and to this end tight optical tolerances 
of one arc-second resolution and less than one arc-second change 
in focus during a one-half hour observation were imposed. It 
was evident that the accurate optical performance analysis of 
tnis sensitive and complex system would require detailed 
structural and thermal models to provide the necessary predic- 
tions of optical component shifts due to thermally induced 
aeforaations. Given that detailed NASTRAN structural models of 
the telescope and the spectrograph already existed, the lUE SI 
Clearly qualified_as an ideal candidate for the first large- 
scale use of the NTA for flight project support.* 


The NASTRAN Structural Model 

TriT:. 2 and 3 are NTA plots of the major structure of the 
lUE telescope and spectrograph (plots of radiating elements were 
the most convenient to include here, but they are essentially 


The NTA was previously used to support the* OSO project 
(reference 8) , but that effort was essentially intended to test 
the integrated analysis concept for debugging purposes. 
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direct overlays of the planar structural elements) . In 
combination these models consisted of 3006 degrees of 
freedom (reference 5 ) , approximately 550 GRID points and 
900 structural elements, and would seemingly be very 
difficult to convert to thermal models due to their sheer 
size and complexity. This problem is specifically 
addressed in the next section, where it will be assumed 
that the reader is familiar with the basic organization of 
a NASTRAN problem, and is aware that the NASTRAN User's 
Manual (reference 3) and the NTA Manuals (reference 1) may 
be consulted for detailed information on the implementa- 
tion of any of the cited techniques. 

GENERAL PURPOSE TECHNIQUE FOR THE CONVERSION OF A 
NASTRAN STRUCTURAL MODEL TO A NTA THERMAL MODEL 

The following Sections are a general purpose state- 
ment outlining the techniques used to convert the lUE 
NASTRAN model to a NTA model. 

Pre-modification 

(1) Acquire a working understanding of the thermal 
design of the structure to be analyzed, with special atten- 
tion being paid to areas with high power dissipations and/ 
or crucial thermal contact couplings. It will often be the 
case that the structural model will be of insufficient 
detail in these areas, requiring additional modelling to be 
done. 


This type of problem is too complex and non-general to 
be discussed in detail in this paper, as it falls into the 
category of engineering judgement; 

(2) Obtain NASTRAN element and GRID point structure 
plots similar to those in Figures 2 and 3, along with a 
verbal description of the model including materials used 
and simplifying assumptions; 

(3) Decide whether NLSS or NLTR results are required 
(linear problems may be solved by the same solution algor- 
ithms, or, for linear steady-state problems, a special 
algorithm is available) . The conversion technique presen- 
ted will initially produce a NLSS model, but the further 
changes required to produce a NLTR model will also be 
described. 
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Card Removal 


Using a punched deck copy of the NASTRAN structural model: 

(1) Remove the Executive and Case Control Decks; 

(2) Remove all Bulk Data Deck cards except GRID, SPOINT 

connection, property, coordinate, relevant $ ,GRDSET BAROR 

S^, BgGIN BULK , and ENDDATA cards (the thermaiiriHipflirl' 
connection cards are listed in the NTA Manual (reference 1) ; 

(3) Remove any permanent constraints defined on GRID cards 

or on a GRDSET card. 

This reduced structural model Bulk Data Deck is now ready 
for conversion to a NTA model. ^ 


NLSS Card Additions 

Using the card deck generated in the previous section: 

(1) Add Case Control and Executive Control Decks; 

(2) If degree-of-freedom information is requested by a card- 

type, always use 1; -a i 


(3) Add MAT4 , MATT4 , MATS , and/or MATTS cards to define 
thermal conductivities (eg - W/m-°C) for each material type 
referenced by the property cards in the model; 

CELASi cards are present, examine them to 
determine, based on their location in the structure, what effec- 
tive conductance (eg-W/°C) should be entered in field 3 of each 
^ a rci ; 


(S) List on an SPCl card all GRID points which 
held at fixed temperatures; 


are to be 


(6) If any GRID points are to be maintained at equal but 
unspecified temperatures, list them on one or more MPC cards 

temperature relationships may be defiJTid if none 
Ot the GRIP points involved are subject to non-linear effects) ; 

IV ^5 Permanent constraints were removed 

from the structural model, carefully examine the GRID points 

-ti.factorlly cSn^Ted into 

(8) Define a temperature "guess vector" consisting of TEMP 
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cards and/or a T EMPO card. All GRID points in the model must 
have a guess temperature/ and those GRID S listed on an SPCl card 
will be constrained to that temperature. In addition, to assure 
convergence, all guess temperatures should be no less than 80% 
of the final steady-state absolute temperature of the respective 
GRID points; 

(9) Define a thermal load vector if heat is to be applied 
to any elements or GRID points; 

(10) If radiative, convective, and/or contact heat exchange 
is to be allowed, CHBDY elements and supporting property and/or 
material cards must be added to define the active surfaces. In 
addition, if gray body radiative couplings are present, RADLST 
and RADMTX cards are required to, respectively, define the 
radiatively active CHBDY cards, and the area-times-view-f actor 
coupling (eg-m^) between the active CHBDY cards; 

(11) Add four PARAM cards (EPSHT, MAXIT, TABS, and SIGMA) to 
define, respectively, the convergence criteria, the maximum 
number of iterations, the absolute temperature conversion factor 
(eg-273.15 for “O , and the Stefan-Boltzmann Constant. 

This completes the basic Bulk Data Deck changes required to 
produce a NLSS NTA deck from a NASTRAN structural model. 


Further NLTR Card Additions 

The NLSS Card deck produced in the previous section may be 
converted to a transient form if the following changes are made; 

(1) Minor modifications to the Executive and Case Control 
Decks are required; 

(2) Proper thermal masses (eg-J/°C) must be applied to the 
GRID points in the model. This is done primarily through speci- 
fying a volume thermal mass (eg-J/m^-°C) on the thermal material 
cards. Additional thermal mass may be applied through the use 
of CVISC and/or CDAMPi cards. The CHBDY card may also be used 
to provide thermal mass in conjuiiction with convection (h of the 
thermal mass specified will be applied to the ambient point(s)). 
Note that structural mass has no relationship to thermal mass, 
with the correct analogy equating viscosity or damping with 
thermal mass. Singular thermal mass matrices are allowed, but 
my experience has been that stability is somewhat impaired; 


*Non-structural mass from the structural model will not supply 
any thermal mass. 
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TL0AD2^card{s^^^"^^^ loads must be referenced via TLOADl and/or 

__ j condition" temperature set composed of TEMP 

cards and/or a T EMPO card must be supplied. The thLmal "qu^ 
vector used in the NLSS model should still be used, as it appears 
matJicertf^^^ stability of problems with singular thermal 

with a aniSS^ oscillating solution, good results are obtained 
with a guess at approximately the average solution value) ; 

(5) are unaffected, but SPCs may only be used to 
constrain GRID points to 0°. They must be replaLd with the 
well known ground-and-load" system* which allows both fixed 
temperatures and fixed time/temperature histories to be specified; 

(6) A TSpP card must be added to specify the number of 

outputr^°'' length, and the frequency of 

vert changes required to con- 

vert a NLSS NTA model to a NLTR one. 


lUE THERMAL PREDICTS vs TEST RESULTS 

These procedures were followed to transform the structural 
m « telescope and spectrograph into a comp^Itr 

9). A comparison at representative 
points of the predicts generated by this model with the actual 
results of an ETU thermal vacuum test (reference 6) performed 
t the Goddard Space Flight Center is presented in Table 1 and 
IS accompanied by the thermal predicts produced by the post-?est 
updated model (reference 10) . Agreement was in gLeral^ScellJ^t 
with the major differences being traced to incorrect coating 

secondary mirror and incorrect 
3 int conduction assumptions at several points in the model. 

ANALYSIS AIDS 


course of using any program extensively, one slowly 
aidS collection of especially effective or uLful analysis 

techniques which through lack of documentation and/or 
lack of general usage are not obvious to the new user. The 
following four items, all used in the lUE analyses, fall into 
this category and should prove useful to most NTA analysts 


*See Reference 1 
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The VIEW Program (references 11 and 12) 

This program is designed to automatically produce the 
raDMTX and RADLST cards required to simulate diffuse radiative 
interchange , including shading, between CHBDY elements. As 
input, it requires_primarily the same CHBDY , PHBDY , and GRID 
cards used in the NTA model, with one additional card type being 
used to specify mesh sizes and optional shading information. 

VIEW is currently operational on IBM 360-370 series computers 
(but could be converted to UNIVAC or CDC operation in roughly 
one month) and executes in as little as IlOK decimal 8-bit bytes 
from an overlayed load module occupying_ 26 7294-byte tracks. 

VIEW was found to be vital to the lUE NTA analyses as it greatly 
simplified the required computation and correction of the large 
view factor matrices involved. 


Modelling Multi-Layer Insulation (MLI) 

A common problem in thermal modelling is how to account for 
MLI. If it is assumed that the heat loss through the MLI cannot 
simply be completely neglected, two methods are suggested: 

(1) Effective emissivity. The NTA assumes that the infor- 
mation entered on the RADMTX cards is area-times-view-f actor 
for gray bodies, rather than Script-AF data. Therefore the use 
of effective emissivity is approximately valid only when the 
CHBDY card using an effective emissivity sees only essentially 
black fixed temperature surfaces (such as the outside of the 
lUE Telescope looking at space) . Also, of course, an effective 
absorbtivity would have to be specified to account for external 
inputs such as solar energy. 

(2) Effective conductance per unit area . This is a more 
straight-forward technique, and requires that the analyst 
estimate the effective conductance through the MLI over a unit 
area (eg-W/m^-®C) and assume that the CHBDY cards are on the 
outer layer of the MLI.* The emissivity and absorbtivity of the 
CHBDY cards would correspond to that of the outer surf ace of the 
MLI, ^and convection would be defined from the same CHBDY cards 
to the GRID points on the structure directly beneath the MLI, 
with the convective coefficient "h" corresponding to the effec- 
tive conductance per unit area through the MLI (see Figure 4) . 
This method is obviously much more flexible than the one pre- 
viously described, as there are no constraints on the surroundings 
nor is there any need to calculate an effective absorbtivity. 
However, care must be taken to ensure that, for both NLSS and 

*New GRID points must be defined for these elements. 
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NLTR runs, the thermal "guess vector" 
grid points associated with the CHBDYs 
ities will result. 


is appropriate for the 
on the MLI, or instabil- 


Solution and Transient Printer Plotting of Linear Equations 

analyses simple, linear parametric equations 
were developed relating the defocus of the telescope to factors 
such as axial and circumferential thermal gradients, changes 
in mirror temperatures, etc. It became apparent that these 
equations could be solved during a NTA run by adding new GRID 

points whose values would be appropriately defined by MPCFl^ It 

would also be possible to use the NTA printer-plotter “if the 
run were a transient, to plot the value of the GRID point which 
represented defocus as a function of time (see Figure 5*) . This 
all'filS^ analysis technique would obviously be applicable to 
linear equations and provides an excellent quick- 

s?ack"Sf 15: “ “"herwlse difficult to peruse 


On-Line NTA Post Processors 

When a large NTA model is to be run more than 10-15 times 
can develop with respect to actually examining the ' 
enormous volume of data which is generated. While transient 

degree, an on-line post processor 
analyses answer, as was found to be the case for the lUE 

quite simple, and involves requesting temper- 
ature, load, and constraint power output in punched as well Is 
printed format (eg-OLOAD (PRINT, PUNCH) = ALL). The punched 

a ^ direct-access storage unit instead of 

a card-puncher, and when the NTA execution is complete, a user 

program may be compiled and executed in-line to 
read the card images and print out neatly formatted summaries 

5°'' 6 shows the format of sample 

pnched temperature, load, and constraint force cards, while 

si^ary report automatically and immediately generated 


*A defocus/time plot from the NTA OSO 
shown as it was the more reproducible 


analyses (reference 8) is 
of the plots available. 
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PROGRAM MODIFICATIONS 


During the course of the lUE thermal analyses, it bec^e 
obvious that several relatively minor improvements to the NTA 
coding would allow a considerable increase in the program's 
efficiency, capability, and "useability . " Since inclusion of 
the actual FORTRAN code changes would be too lengthy for this 
presentation, discussion will be_limited to an outline of the 
use of the new features and the NTA subroutines which were modi- 
fied to implement them, with the understanding that more complete 
information can be obtained from either Mr. Clifton E. Jackson, Jr. 
or Mr. Reginald S. Mitchell, Code 720.2, Goddard Space Flight 
Center, Greenbelt, Maryland 20771. 


Additional "Space Node" Flexibility 

An option was supplied to allow elimination of the automatic 
assumption that any energy not accounted for in the RADMTX should 
be discarded (the so-called "space node" feature) . The new 
option, if elected, would assume that energy not accounted for 
would be returned to the element from which it was radiated. 
Additional related features added included the capability to 
partition the RADMTX into "space node" and "non-space-node" 
portions, and a summation check to assure that the RADMTX does 
not in effect specify emissivities greater than 1.0. The 
addition of these features required modifications in subroutine 
RMG only. 


Script-AF Input 

The original’ NTA required the use of an internal Scrijpt-AF 
generator in order to include radiative interchange, a fadt that 
prevented the use of externally generated Script-AF matrices. A 
modification to subroutine RMG now optionally allows either 
symmetric or unsymmetric Script-AF matrices to be input to the 
NTA using the standard RADMTX and RADLST cards. 


Symmetric Decomposition 

When the full-scale NLSS model of the lUE experiment was 
first run, CPU and I/O solution times on the order of an hour 
were observed. An examination of the run logs for this problem 
indicated that large amounts of time were being spent in the 
unsymmetric decomposition of symmetric radiation and conductance 
matrices. It was determined that minor limitations on input 
(emissivities not equal to exactly 1.0, and uniform thermal guess 
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vectors for unconstrained GRID points if radiative interchange 
were included) would allow the use of symmetric decomposition 
in both cases. Modifications were made to subroutines RMG and 
SSGHT, an additional parameter was added, and run times were cut 
to less than 18 minutes CPU and I/O, with only round-off error 
changes in the answers produced. A similar conductance matrix 
decomposition feature was later added for transient runs, and it 
appears that it would be especially useful if time step sizes 
were changing frequently. 


CONCLUDING REMARKS 

The NTA has been used to provide thermal analyses of the 
lUE SI using structural models as the analytical baseline. The 
verification of this technique via test has demonstrated the 
accuracy of the NTA as a thermal analytical tool, and the analysis 
aids and program modifications listed will assist in making it 
a useful one. 
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TABLEI - CORRELATION OF THE lUE TEST AND THE ANALYTICAL MODELS 
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Figure 2. NTA Plot of the Major lUE 
Telescope Radiating Elements 
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Figure 3. NTA Plot of the Major lUE 
Spectrograph Deck’s Radiating Elements 
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Figure 5. NTA Plot of Defocus (vs) Time for the OSO Telescope 
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Figure 6. Partial lUE Punched Temperature, Thermal Load, and Constraint Force Output 
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Figure 7. Sample Automatic Post-Processor Summary Report for the lUE SI 


